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Abstract. We analyze electron spin relaxation in electronic transport through 
coherently coupled double quantum dots in the spin blockade regime. In particular, 
we focus on hyperfine interaction as the spin relaxation mechanism. We pay special 
attention to the effect of the dynamical nuclear spin polarization induced by the 
electronic current on the nuclear environment. We discuss the behaviour of the 
electronic current and the induced nuclear spin polarization versus an external 
magnetic field for different hyperfine coupling intensities and interdot tunnelling 
strengths. We take into account, for each magnetic field, all hyperfine mediated spin 
relaxation processes coming from the different opposite spin levels approaches. We 
find that the current as a function of the external magnetic field shows a peak or a dip, 
and that the transition from a current dip to a current peak behaviour is obtained by 
decreasing the hyperfine coupling or by increasing the interdot tunnelling strength. We 
give a physical picture in terms of the interplay between the electrons tunnelling out 
of the double quantum dot and the spin flip processes due to the nuclear environment. 
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1. Introduction 

In the last decade solid-state spintronics and quantum computing have experienced a 
great development. In particular, quantum dots (solid-state fabricated zero dimensional 
devices) have been widely investigated both experimentally and theoretically. For 
quantum computing and quantum information they have become major candidates for 
implementing quantum bit units [1] (qubits), but also from a fundamental point of view, 
since quantum dots resemble artificial atoms, they are highly interesting systems for 
studying basic atomic physics. In this context, spin decoherence and relaxation are 
among the most desirable mechanisms to be understood since they represent the main 
sources of quantum computing errors. 

We investigate spin relaxation in double quantum dots (DQD). Spin Blockade 
(SB) [2j is a very suitable regime for attempting single electron spin manipulation 
because two electrons are trapped in a DQD, since Pauli Exclusion Principle avoids 
electron transport through the dots. Moreover, SB is attainable through transport 
experiments in DQDs in several materials. However, spin relaxation processes can 
partially destroy the SB releasing the trapped electrons, leading to a small, though still 
measurable leakage current (~ pA-fA) [2HB]. Spin-orbit coupling [7], cotunnelling [HIE] 
and hyperfine (HF) interaction between the DQD electrons and the surrounding 
nuclei spins of the host material represent the main mechanisms for spin-relaxation. 
Depending on the material, one or more mechanisms may be involved collaborating or 
competing [3HUE]. 

In this paper, we study theoretically spin relaxation in a DQD in SB regime due 
to HF interaction with the lattice nuclei spins, which has been in the last years a very 
active field both experimentally [31H1I61 IT014T7] and theoretically [T8H23] . We pay special 
attention to the dynamical nuclear spin polarization induced by the electronic leakage 
current [T3l[T8l 12211211 125] emerging from the spin-relaxation transitions (figured]). This 
effect is often not taken into account when studying the current through the DQD in 
SB regime. 

In the present work we calculate both the electronic leakage current and the nuclei 
spin polarization induced by the electrons tunnelling through the DQD. In the SB 
regime, two electrons in the DQD can be either in a triplet or in a singlet state. 
Nevertheless, current is only allowed to pass through the DQD when they are in a 
singlet, otherwise they remain trapped in a triplet state. HF interaction mixes triplet 
and singlet subspaces, and thereby lifts SB. The mixing is due to the different HF 
interaction strengths, i.e., two different effective magnetic fields (induced Overhauser 
fields on the electrons by the nuclei), within each dot. In addition, scattering processes 
between electron and nuclei spins that lead to spin relaxation, also induce a non- 
negligible nuclear spin polarization as the current flows through the DQD. Moreover, 
this induced nuclear spin polarization itself, as we shall see below, modifies the mixing 
between the triplet and singlet subspaces, acting back on the electronic current through 
the DQD. Furthermore, we find that the Overhauser field proportional to the electron 
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(a) DQD, leads and the lattice nuclei. 



(b) Transport scheme. 



Figure 1: (Colour online.) a) A coherently coupled DQD is coupled through tunnelling 
barriers to leads and to the surrounding nuclei of the host material, b) Transport 
scheme. Electrons tunnel from the left lead into the DQD states (solid blue arrows 
online). When electrons fall into the singlet states (S±), tunnelling out of the DQD to 
the right lead is allowed (solid blue arrows online). However, when electrons fall into 
a triplet state (T±) Pauli Exclusion principle prevents them from tunnelling out of the 
DQD into the right lead, thus, SB occurs. Finally, if electrons tunnel into the T x state 
(a mixture of a singlet state and the antiparallel spins triplet T ), they can only tunnel 
out of the DQD to the right lead if there is a net nuclear spin polarization (dashed blue 
arrow online), otherwise T x becomes T , which is also a blocked state. Electrons trapped 
in T± states interact with the nuclei and relax to a singlet state or to the T x state (solid 
red arrows online) through spin-flip processes. However, electrons in T x state can also 
spin-flip back to a SB state (red dashed arrows online). Therefore, when electrons are in 
the T x state there is a competition between the tunnelling rate from the T x to the right 
lead, and the spin-flip rate from T x to the T± triplets (a detailed discussion is given in 
section [2]). We will show that this competition will give rise to different physical features 
in the tunnelling current as a function of an external magnetic field (see section [3]). 



current induced nuclear spin polarization, is in general larger than the one obtained 
when just considering HF interaction as a random stationary magnetic field acting on 
the electron spins. 

In previous works [2211211 EH], the interdot tunnel coupling between the quantum 
dots was considered incoherent, namely, the system was assumed to be in the sequential 
interdot tunnelling regime. However, in the present work we focus on the resonant 
tunnel regime and we consider coherently coupled quantum dots (figure fTal) . We have 
found in this case a different behaviour both for the current and the induced nuclear 
spin polarization with respect to previous works. In the coherent coupling regime, in 
addition to the two spin parallel triplets, there is a spin antiparallel triplet. Thus, 
in addition to the spin-flip transitions between the singlet and triplet states, spin-flip 
processes between triplet states are also present, and contribute to the nuclear spin 
polarization and the electronic leakage current. Furthermore, in the present model we 
find that the transition rates between the DQD and the leads depends on the nuclear 
spin polarization (figure fib]) . This effect will lead to new features in the current through 
the DQD. 

The paper is organized as follows. The Hamiltonian for the DQD coupled to the 
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leads and the surrounding nuclei spins, and the rate equations for the occupation of 
DQD levels and the nuclei spin polarization, are considered in section |2j The results 
and discussions are considered in section [3j And finally, the conclusions are presented 
in section HI 

2. Model 

2.1. Hamiltonian and DQD eigenstates 

The system we investigate is a DQD coupled to two uncorrelated electron reservoirs 
(leads) and to the nuclei spins of the surrounding host material (figure [Tall . We consider 
a spin up and a spin down level in each dot. The Hamiltonian is the following: 

H = H DQD + H leads + V LR + V T + V HF (1) 

where Hdqd and Hi ead s are the Hamiltonians for the isolated DQD and the electron 
reservoirs, respectively, Vlr is the inter-dot tunnelling Hamiltonian, and Vp and Vhf 
correspond to the DQD coupling with the leads, and the HF interaction with the nuclei 
spins, respectively. We neglect cotunnelling processes because we consider that the 
tunnelling coupling through the contact barriers is much smaller than the thermal energy 
and the bias voltage [26]. Moreover, the energy of the DQD levels are strongly detunned 
respect to the contacts chemical potentials [HIE] - 

We consider contact HF interaction, as we regard electronic wave functions with s 
like symmetry [27]. Hdqd and Vlr are: 

H D qd = e i™i° + U i fl n fl H + u lr ^2 + 9^BB cxt Si z 

la I aa' I 



V L r = t L R ^2 (d-La^Ra + h.c 



where I = L (left dot), R (right dot) and a =t, 4- d\ a {di a ) creates (annihilates) an 
electron with spin a and energy e\ in the Z-th dot. hi a = d\ a di a is the occupation 
number operator and Si the electron spin operator. U\ {Ulr) is the intra-dot (inter-dot) 
Coulomb interaction, 5 ext the external magnetic field and tLR the tunnelling matrix 
element between the dots. We do not focus on any particular material, thus we take 
g = 2. H i eads and V T are: 

Hleads = ^ ^-Ika^lhc/^lka 
Ika 

V T = J2(llk4kadla + h.C.) (2) 
Ika 

where c lka . (cn- a ) creates (annihilates) an electron in the Z-th lead with momentum k, 
spin a and energy e ika . 7^ are the tunnelling matrix elements between the dots and the 
contacts. Finally, the HF interaction term Vhf is: 

N 

(3) 

l=L,R i=l 
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where A\ = zM|\I/o(rj)| 2 is the HF coupling [23] between the i-th nuclei spin I f at site 
rj and the electron spin Si, where v is the volume of a unit cell containing one nuclear 
spin, A characterizes the hyperfine coupling strength and ~§>\{ji) is the single-particle 
electronic wave function for dot /, evaluated at site r;. For simplicity, we only consider 
spin- 1 nuclear species. The HF couplings for the left and right dots depend on the 
square modulus of the electron wave functions at the position of the nuclei, therefore, 
for realistic dots they are in principle different for each dot [29||30]. We consider the 
case of homogeneous hyperfine couplings [31J within each dot, hence, A\ = Ai/N, where 
iV is total number of nuclear spins, and split Vhf (see (jlj)) in the following two terms: 

Vhf = J2i=l,r # Eli (l (SiJi- + Si-h+) + S lz I lz ) = V sf + V z , (4) 

where S± = S x ± iS y and I± = I x ± il y are the raising and lowering operators for 
the electron and the nuclei spins, respectively. V s f corresponds to the x-y-components, 
which are perpendicular to the external magnetic field and, thus, is responsible for the 
flip-flop transitions between electron and nuclei spins; and V z corresponds to the z- 
component, which is parallel to the external field and, hence, contributes to the Zeeman 
splitting, as we will see below. We treat these two terms separately. For V z we perform 
a mean-field approximation which gives |22j : 

v z MF = IE^Ap (5) 

where P = ((N^) — (N±))/N is the net nuclear spin polarization, and (iV-j-) ((A^)) is 
the average number of spin up (down) nuclei. V^ 1F is now equivalent to an effective 
magnetic field within each dot induced by the nuclei on the electrons (the Overhauser 
field), and proportional to the nuclear polarization given by: 

B L(R) = A L{R)P ( , 

This effective field is in general different for each dot |19j . and gives rise to an effective 
Zeeman splitting within each dot that adds to the one produced by the external field. 
Below we shall see that this inhomogeneity is essential for lifting SB [3]. We consider 
Vt and V s f as perturbations. Thus, the unperturbed Hamiltonian for the isolated DQD 
is the following: 

= Hdqd + Vlr + V Z MF . (7) 

SB occurs when the source-drain voltage is tuned so that the number of electrons in 
the DQD varies between one and two. The right quantum dot is always occupied with 
one electron while another electron can tunnel from the source to the drain through the 
DQD. The transport scheme is the following: 



(0,a) — ► (a» 




o = o' 



where in (n, m), n (m) accounts for the population of the left (right) dot level, and 
a, a', a" =t, 4- Interdot tunnelling (Vlr) is allowed only between states with the same 
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total spin. Since the state (0,44) has total spin zero, while the states with = 0' have 
total spin one, the transitions (a, a) — > (0, 44) are forbidden. The Hilbert space that 
we have considered, consists of the atomic basis {|0,t), |0,4), | 4,4), I 4>t)> I t>t)> 
| 4,4), 1 0,44)}- We investigate the coherent resonant transport regime. In this regime, 
the energies of the DQD two electrons states \a, a') and |0, 44) are degenerate in the 
absence of a magnetic field, i.e., the so called zero detuning regime. In this case Hq 
(see ([7])) is exactly diagonalizable. Its eigenenergies and eigenstates are the following: 

E Tx = e L + e R + U LR 

E s± =e L + e R + U LR ± V2 Mt LR 



E T± =e L + e R + U LR ± (^gn B B c ^ t + 4p 



(9) 



and: 



where: 



|T + ) = |t,t> 
|T_) = |4,4) 
\T X ) = ± (\T ) - x\S 02 )) 

\S ± ) = -L f\ Sn ) ± 1 (\S 02 ) + rr|T ») (10) 
|T ) = ^(|4,4) + I4,4)) 

|5ii) = ^ (I 4,4) -I 4,4)) 

|5 02 ) = |o,44) (11) 



being: 



N = VTT 



X 



2 



1 A_P 



x 



v/2 2t LR 



A ± = ^(A L ±A R ). (12) 

Thus, the two electron molecular basis is {|±), \S±), \T X ), \T±)}, where | + (— )) = |0,4 
(4)) are the single electron states. Equation (TTUT) shows that Vl R mixes Su and S02 
singlets, and V^ IF mixes T triplet with Su and 602 singlets when A L ^ A R . The singlet- 
triplet (ST) mixing is given by the weight xjM . This quantity is the ratio between the 
Zeeman splitting difference within each dot (A_P/2), and the exchange energy defined 
as \Eg ± — EtJ = \/2MtL R . Therefore, the ST mixing depends on the competition 
between these two energy scales. A large (small) difference between the HF coupling 
intensities, and small (large) interdot tunnelling strength increases (decreases) the ST 
mixing (see ffTU]) and (fT2"l) ). Furthermore, x is zero when either the HF couplings have 
the same value for both dots (A L = A R ) or the nuclei spins are completely depolarized 
(P = 0). In both cases the usual singlet-triplet basis is recovered. Due to the mixing 
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with T triplet state, S± are not pure singlet states anymore, however, for simplicity we 
shall continue calling them singlet states. Finally, notice that now the eigenenergies of 
singlet states also depend on the nuclear spin polarization through M (see (jHJ) and (|T2l)). 
The transport scheme in the molecular basis is the following: 



This transport scheme shows that for the coherent inter dot tunnelling regime, when 
having different Overhauser fields in the dots, a current channel ((0, a) — > T x — > (0, a')) 
is opened. Therefore, there are two possible situations: i) x = 0. In this case the 
incoming electron will fall either in a singlet state {S±) or in a triplet state (T or 
T±). Thus, there are two transport channels and three SB states; ii) x ^ 0. The 
incoming electron will fall either in a singlet state (S±), the T x state or in the T± 
triplets. Therefore, now there are three transport channels and two SB states. However, 
once the electrons drop in a SB state, the current drops to zero. 

2.2. Rate equations 

The Hamiltonian (CQ) can be written now as follows: 



between the leads and the DQD, namely, between one-electron and two-electron states. 
V s f is responsible for the spin flip-flop transitions between the DQD electron spins and 
the surrounding nuclei spins. The time evolution of the DQD molecular states are 
obtained with the following rate equations: 

pt± = w T±)S+ ps+ + w T±:S _ps- + w T± , Tx p Tx + r T±i± p ± 

-(W s+tT± + W S „, T± + W Tx , T± )p T± 

p Tx = W Tx ,T + pT+ + W Tx ,T.pT. + T Tx ,+p+ + T Tx -p- 

-(w T+ , Tx + w T _ iTx + r+, T!C + F-,T X )PT X 
ps ± = w s± ,t+pt+ + w s± ,t-pt_ + r s± ,+p+ + r s± -P- 

-{w T+ , s± + w T _, s± + r+, 5± + r_ iS± )p 5± 
p± = ^±,s+ps+ + ^±,s-ps- + f±,t x pt x - (r s+ ,± + rv i± + r Tx ,± + r T±i± )p ± (15) 

where is the occupation of the i-th state. Tij (Wj,/) are the tunnelling (spin-flip) rates 
between an initial DQD state \i) and a final DQD molecular state |/). Both tunnelling 
and spin-flip rates are computed in section 12.31 

The spin electron-nuclei flip-flop processes induce a non-negligible nuclear spin 
polarization, which will be positive or negative depending on the specific spin-flip 




(13) 
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processes. In the following scheme we show the spin-flip processes which contribute 
to each nuclear spin polarization direction: 

T + ^{S + ,S„,T x }^P>0 
T..^{S + ,S-,T x }=>P<0 
{S + ,S^,T x }^T + ^P<0 

{S + ,S-,T x }^T_=>P>0. (16) 

E.g., the process T + — > {S + , S-,T X } flips down an electron spin and up a nuclear spin, 
thus, polarizes positively the nuclei spins. Furthermore, the nuclear spin polarization 
becomes dynamical due to the electrons tunnelling through the DQD. Therefore, we 
describe the time evolution for the induced nuclear spin polarization using the following 
rate equation: 

P = (W T _,s+ - W T+ ,s+) Ps+ + (W T _,s. - W T+ ,s.) PS- + (W t _,t x - W T+ , Tx ) Pt x 

+ {W S+ ,T+ + W S _,T+ + W Tx ,T + )pT+ ~ {W S+ ,T- + W S _,T. + W Tx , T _)p T _ 

-W vcl P (17) 

where W re \ is a phenomenological rate that accounts for the nuclear dipole-dipole spin 
interaction which is responsible for nuclear spin depolarization. In this equation we have 
assumed that in absence of spin-flip processes (W/,i all zero), the nuclei spins completely 
depolarize, namely, that temperature is much larger than the nuclei spin level splittings. 
As will be shown in section 12.31 both the tunnelling and the spin-flip rates depend on 
the nuclear spin polarization. Therefore, ( JT5l) and ( JT71) form a set of eight non-linear 
equations which we must be solved numerically (section [3]). 

2.2.1. Reduced model In order to get physical insight on the results which will be 
discussed in section [3j we have developed a simplified model for the rate equations ( |T5l) 
and (TTTj) through the following assumptions: i) we consider that spin-flip processes are 
effective when the electrons are in a SB state. Thus, once the electrons are in a singlet 
state it is much more probable for them to tunnel out of the DQD than to flip its spin 
to a T± triplet state. Therefore, we neglect the singlet to triplet spin-flip transitions 
(Wt±,s+ an d Wt±,S-)- h) The singlet states and the one electron states empty much 
faster than the T± triplets (SB) and the T x state. Thus, in the time scale on which spin- 
flip transitions are relevant we assume that the occupation of the one electron state and 
of the singlet states have reached their stationary value, namely, p± ~ ps ± ~ 0. iii) 
The nuclei relaxation time due to the spin dipole-dipole interaction is long enough to 
be neglected (W Te \ —> 0). Under these conditions, the rate equation for the nuclei spin 
polarization ( IT7|) is related to the rate equations for the triplet states (fT5|) through the 
relation —pT + = P, thus, pr_ —pr + = P (where we consider pt_ (0)—pT + (0)— P(0) = 
and p±(0) = ps ± (0) = as initial condition). Taking these considerations into account 
the rate equations (115]) become: 

p T = (2u% + 0) p Tx - « + JDpt + {oj s 1 + u*)P (18a) 
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p Tx = « + u'Dpt - {u s l + u;*)_P - (2w* + 0)p Tx (18b) 
P = 2u*pr m + (wf? + w")p T - «' + oj%)P (18c) 



where p^ = px + + Pt„ and 
1 



1 (^5 + ,t + + Ws.,r + ) ± (ty s+ , r - + W s _, r _) 

2 i + r/r +>5+ 



or 

We have considered the same value of the tunnelling coupling for both contact barriers, 
T L = T R = T (see section I2.3.2p . Finally, by summing ( I18ap and fj!8bj) we find that 



Pr + Pt x = 1, where we consider as initial condition: pr(0) + Pt^(O) = 1, thus, the rate 
equations become: 

p Tx = (w* + <) - (3w* + /? + Or, - (w* + 

P = (w« + + (w tt - u;f)p T , - (w* + uf)P (20) 

The rates u± account for the SB lifting due to ST spin relaxation: T± — > {S+, £_}. u± 
((3) account for the spin-flip (tunnelling) rates from T x to T± (T x to |±)). In section 12.3.21 
we will see that the rate ^+,s + — r/2, hence, it can be regarded as constant. In section[3] 
it will be shown that the current through the DQD depends on ratio between and 
0. The following scheme offers a qualitative description of this dependence for the case 
of A_ ^ 0: 

— y (0, a) — > High Current Regime 
U +>P T± — > Low Current Regime. 

When u;" ^> pi electrons in T x state have a larger probability to spin-flip to the T± 
triplets than to tunnel out of the DQD, and only the leakage current coming from the 
T± — > {S + , S*_} transitions provide a leakage current, which we define as the low current 
regime. On the contrary, when w+ <C (3 electrons in T x state have a larger probability 
to tunnel out of the DQD than to spin-flip to the blocked T± triplets, and the current 
is enhanced. We define this the high current regime. In section [3j we will see 

that the competition between spin-flip and tunnelling transition rates determine the 
behaviour of the leakage current. Moreover, we will see that the transition between 
the low current regime and the high current regime is obtained by varying either the 
interdot tunnelling strength or the HF coupling intensity. 



2.3. Transition rates 

2.3.1. Spin- flip rates The transition rates between H eigenstates (see (JTOl ) due to 
spin-flip processes are calculated in perturbation theory by means of Fermi's Golden 
Rule [20ll22p26j . The left (right) dot electron spin z-projection is given by m,L{R) = ±1/2, 
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hence, the total spin projection in the ^-direction for a DQD two electron state is 
M = mi + mn (M = —1,0,1 for triplets, and M = for singlets). The spin-flip 
interaction term V s f (see (J4|)) increases an electron spin by one while decreasing a nuclei 

y g j, y g j, 

spin by one also (and vice versa), thus, M : —1 1. We consider the initial 
state |ijv)|aA/) which consists of the initial nuclei states \i^) and the DQD electron states 



\a M )- \in) is given by \i N ) = \m h m 2 , 



,rrij, 



,m N ) 



, where nij 



±1/2 is the spin 



of the j-th nuclei, and G {|5±), l^), |T±)}. The final state \ fn)\^M') is connected 
to the initial one by having the j-th nuclear spin flipped and a different electronic state 
|/?m') with \M — M'\ = 1. Therefore the spin-flip rate for the transition that flips up an 
electron spin and down a nuclei spin is: 

N 



w, 



27r^^|(/3 M+1 |(/ JV |V; / |i JV )| aM )| 2 >V ijv 5 (E Pm+1 - 5, - E aM ) 



(22) 



2tt 



1 1 



Ai([3 M +i\Sf\a M ) 



l=L,R 



N 



^N\^Ii\i N ) m N 5 {E Pm+1 - 5, - E aM ) 

3=1 *iv 



where |/jv) = jj"|ijv) and we have taken H = 1. The sum over initial states runs over 
all configurations of the internal degrees of freedom, i N , that give the state \in). Each 
state is weighted by the probability of having that configuration, which is given by the 
distribution function Wi N . E aM (Ep M+1 ) is the energy of the initial (final) electronic 
state \aj[,f) (\Pm+i)), an d Sj is the energy splitting between the up-down levels of the j-th 
nuclei spin. We assume independent nuclei spins, hence, Wi N = W mi x VV m2 x . . . x W mjv , 
where W mj is the probability that nuclear spin j has the value rrij. Thus, the sum over 
initial nuclear states becomes: 



mj\Ip] 



\rrij) W„ 



in 



-1/2 



(23) 



1 to eliminate 



where we have used the normalization condition W mfc =i/2 + VV mfc = 
all other sums over mk except rrij. The probability of having nuclear spin j in a certain 
state is related to the overall nuclear spin polarization as: 

(iV» _ 1 + P 



=1/2 



-1/2 



N 



i-p 



(24) 



N 2 

In general, the g- factor is much smaller for nuclei spins than for electrons [32J. Therefore, 
under experimental conditions [2| |4TI6| [12] . the nuclear splitting is usually negligible 
compared to temperature and the energy difference between electronic levels, hence, Sj 
V j — 1 . . . N can be safely neglected. Thus, the spin-flip rate becomes: 

2 

M0M+i\S?\a M ) 



w. 



l=L,R 



5 \ E 0U+l ~ E *m) ■ 



(25) 
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Repeating the same procedure we have that the spin-flip rate for the transition that 
flips down an electron spin and up a nuclei spin is the following: 

2 

l —^5{E PM _ 1 -E aM ). (26) 



7T 



^ A^pM-ilSi \a M ) 



l=L,R 

Notice that these spin flip-rates depend on how much the nuclei are polarized. When the 
nuclei spins are fully polarized in the positive (negative) direction Wp M _ iaM (Wp M+iaM ) 
vanishes. 

We see that the derived spin-flip rate requires energy conservation, so strictly 
speaking leads to zero spin-flip for E aM ^ Ep M±l . However, in reality it is possible 
to exchange energy [33] with the environment e.g. as phonons. We model this by 
replacing the Dirac delta by the following expression: 

1 7 

5 ( e Pm±i ~ e <*m) -> -7— — ~2 r x Cp M±1 a M (27) 

n (Ep M±1 -E aM ) + 7 2 

where: 

J 1 if Ep M±1 > E aM ->■ Energy emission 

Cfai±lQM= \ exp( ^r ga " ) if E Pm ±1 <E* m Energy absorption (2S) 

where T is the temperature and ks the Boltzmann constant. The Lorentzian is maximal 
for the elastic case and falls off with increasing energy exchange on the characteristic 
scale 7. This parameter is assumed to be of the order of the typical phonon energy, 
7 ~ yueV [33] . The function Cp M±l0lM accounts for the low temperature energy 
emission/absorption asymmetry [33 s 34], i.e., it is much easier to emit than absorb energy 
from e.g. a phonon bath. Formally, one can include the electron-phonon coupling as 
a perturbation together with V s f and thereafter use a T-matrix approach to find the 
phonon-mediated HF spin-flip rate [35H37] . However, here we do not pursue an exact 
modeling of the way the energy is exchanged with the environment, but simply include 
the fact that the spin-flip rate decreases as the energy involved in the flip-flop processes 
increases [20,122] • Therefore, the spin-flip rates we obtain are the following: 

2 

1±P 7 



2J A((3M±i\S^\a M ) 



l=L,R 



, x2 „ C /?M±ia M - (29) 

\ E p M ±x- E * M ) +7 2 

In section [3], we will show that the amount of induced nuclear spin polarization depends 
on the competition between emission and absorption processes. The energies that appear 
in the function Cp M±lD , M (see fl28l) ) are the eigenenergies shown in (jUJ). Thus, the energy 
differences in Cp M±1(XM depend on the interdot tunnel strength, the HF coupling and the 
external magnetic field. Therefore, for a given temperature the emission/absorption 
asymmetry can be controlled through these parameters. Finally, the only matrix 
elements between the DQD states which are different from zero in (|29l) are: 

(T x \ (A L St + A R S±) \T T ) = ^± (30a) 
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(S ± \ (a l S~ l + A R S R ) \T + ) = - (a_ T (30b) 

(S ± \ (A L S + L + A R S+) |T_) = A_ ± ^± (30c) 

From these expressions we distinguish between two different HF mediated spin relaxation 
processes: i) the triplet-triplet relaxation ( I30aj) . and ii) the singlet-triplet relaxation 
( f|30bl) and (I30cj) ) . Notice that if the HF coupling intensities have the same value 
for each dot (A_ = 0, and then x — 0) only the triplet-triplet spin relaxation survives 
(T T — > T x = T ) while ST relaxation probabilities become zero. Thus, SB lifting involves 
having A_ ^ 0, as discussed in section [27T1 (see (fl3j) ). Finally, the matrix elements also 
depend on the interdot tunnel through the ST mixing parameter x (see ffT21) ). 

2.3.2. Tunnelling rates The tunnelling rates between the leads and the DQD are 
calculated using Fermi's Golden Rule [261 El]- Given an initial DQD state with n 
electrons \a n ) the tunnelling rate for an incoming electron to the final DQD state |/3 n +i) 
with n + 1 electrons is: 

r^ +1 ,a n = rj(/i D -/iO^|(/Wil4>n>r. (31) 

(7 

And the tunnelling rate for an outcoming electron to the final DQD state \f} n -i) with 
n — 1 electrons is: 

=T l (l-f( t i D -^))J2\Wn~i\dM\ 2 (32) 

where I = L,R. Ti = 27r |7/ fc | 2 I);, where it is assumed that the density of states in 
both leads D\ and the the tunnelling couplings 7^ (see (T5])) are energy- independent. 
f(fiD — Hl(r)) is the Fermi distribution function for the left (right) lead, Hl(r) is the 
chemical potential of the left (right) lead and \1b is the DQD chemical potential. The 
DQD states that appear in the matrix elements in the tunnelling rates are the eigenstates 
of H (see f lTUj) ). The chemical potentials (/i/) are tuned such that the system is in the 
SB regime, and we define Tp n±ltan = J2i=l r ^„ ±1)tt „' Therefore, the tunnelling rates 
different from zero are: 



rv ± ,± 


= r L 








r s+>± 


r L 
4 


(■* 


x 2 

Jp 


2x 

t m 


rv,± 


_ r L 
4 




x 2 


2x 
± M 






1 







LT - ± = YW (33) 

for electrons tunnelling from the leads into the DQD, and: 
r - r - Vr 1 

r ±A = T KJT2 ( 34 ) 
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for electrons tunnelling out of the DQD to the leads. As pointed out in section \2.2\ 
the tunnelling rates also depend on the nuclear spin polarization. Additionally, T± T a 
becomes zero if x = (see (1121) ). i.e., when the HF couplings within each dot are the 
same (A- = 0), or the nuclei spins are completely depolarized (P = 0). 



where I = L,R. 

3. Results and discussion 

In this section, we present the results obtained by solving numerically the full system 
of rate equations derived in section |2] ( (TI5]) and (TTTfl ). The system has three main 
relaxation time scales: i) the tunnelling through the contact barriers, ii) the electron- 
nuclei spin- flip, and iii) the nuclear spin relaxation (oc W re i). The nuclear relaxation 
time is the largest one (~ minutes [3JE1EE])- We are mainly interested in investigating 
the system on a time scale for which the tunnelling and the spin-flip are the relevant 
relaxation times, whereas the nuclear spin relaxation time is much larger. 

We will discuss the behaviour of the electronic current and the induced nuclei 
spin polarization versus the external magnetic field for different HF couplings and 
interdot tunnelling strengths. The HF interaction constant is not well known for all 
materials, thus, we have considered values for the HF coupling intensity in the range 
Al = 70 — 90 fieV [HEIEH]. The difference between the HF couplings in each dot is 
held constant for all cases (A R = 0.8A L ). The plots described below have been obtained 
by solving numerically the rate equations derived in section [2j ( IT5|) and ( |T71) . sweeping 
the external magnetic field from negative to positive values. We have performed the 
calculations for three cases: i) In order to observe only the effect of the dynamical 
nuclear spin polarization, we consider that the nuclei spins are initially fully depolarized 
for each value of the external magnetic field. Thus, the sweeping rate is much slower 
than the nuclear spin relaxation rate W Ie \ (see ( jTTl) ). Nevertheless, for each value of 
the external field the rate equations are solved in a time scale much smaller than the 
nuclear spin relaxation rate in order to capture only the tunnelling and spin-flip events, 
ii) Initially the nuclei spins are completely depolarized, however, while the magnetic 
field is swept, nuclear spin polarization is built up. Thus, the sweeping rate is much 
faster than the nuclear spin relaxation rate, and a feedback processes between nuclear 
spins and electronic current occurs, iii) we proceed as in ii) sweeping the external field 
forwards and backwards. In this case we observe hysteresis, as has been widely observed 
experimentally for different DQD devices. [3| 1^16} ITU] 

As stated in section |2l we consider the zero detuning case (i.e., \a, a') and \So 2 ) 
are degenerate). In previous works [SUES], finite detuning was considered, therefore 
only transitions close to one ST level crossing are important. Now, different energy 
levels approaches participate and we take into account all HF mediated electron spin 




(35) 



n a n ,/3 n ±i 
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Figure 2: (Colour online.) Schematic energy level arrangement for the DQD eigenstates 
in the different external magnetic field regions that are investigated, a) B ext < Bs ± ,t t , 
b) B s±> t t < B cxt < B TT , c) B TT < B cxt < B s±iT± and d) B S± , T± < B ext . 



relaxation processes, at a fixed external magnetic field. Figure [2] shows the energy 
levels for the different regimes considered. -Bs ±i t t (Bs±,t±) corresponds to the value of 
the external field for which the S±T T (S±T±) crossings occur, and Btt the value for 
which the T±T X crossing occurs. Notice that the different S±T T crossings occur for the 
same value of the external field (see Q), and the same happens for the different S±T± 
crossings. The energy emission spin-flip rates which contribute for each range of the 
external field in figure [2] are given in the following scheme: 

B cxt <B s±:TT W s± , T .,W TxtT _,W T+jTx P<0 (36a) 

Brr < < B S± , T± { Ws _^J^ Z III ^ 

B s±tT± <B ext => W s± ,t + ,W Tx , t+ ,W t _, Tx => P>0 (36d) 

where the sign of the nuclear spin polarization in each case is described in ( TT6l) . Only 
emission spin-flip rates are shown because at low temperatures they dominate over the 
absorption rates, hence, they will determine the sign of nuclear spin polarization. Thus, 
depending on the intensity of the external magnetic field there will be different spin-flip 
processes dominating. 



3.1. Dependence on the hyperfine interaction intensity 

3.1.1. Small magnetic fields In this section we consider small external magnetic field 
intensities. We will focus on the Bs ± ,t^ < B cxt < Bs ± ,t ± range of magnetic fields. The 
main spin-flip transition rates corresponding to this range of magnetic fields are shown 
in ( I36bl) and ( 136cl) , and the energy levels scheme is plotted in figure [2b] and figure [2c] 
(see also figure HJ). In this case, the transitions around the T±T X crossing give the 
main contribution to the current, although other possible opposite spin level transitions 
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participate (see (136bj) and (j36cj) ). The S±T T and S±T± crossings occur for larger 
magnetic field intensities and will be studied in section l3.1.3l Figure [3] shows the induced 
nuclear spin polarization and the leakage current through the DQD versus the external 
magnetic field for different values of the HF coupling (At), in this range of the external 
field. 

The induced nuclear spin polarization versus the external field is shown in figure l3al 
(without feedback), figure [3b] (with feedback) and figure [3c] (sweeping forwards and 
backwards with feedback). Figure H] shows the energy levels. When B ext < Btt 
(B ext > B TT ) the net induced nuclear spin polarization is negative (positive): Emission 
processes are stronger than absorption processes, hence, they mainly determine the 
sign of the nuclear spin polarization. In f!36bj) it is shown that when B cxt < Btt 
(B ext > Btt) the dominating spin-flip rates are Ws_,t_, Wt x ,t_ and Wt + ,t x {Ws±,t + , 
Wt x ,t + and Wt_,t x ), which polarize the nuclei negatively (positively). Nevertheless, 
in this regime absorption processes have the effect of partially compensating emission 
processes resulting in a finite but not complete spin polarization of the nuclei. In 
figure [3a] figure [3b] and figure [3c] it is shown that increasing the HF coupling increases 
the nuclei spin polarization. In this region, since \Es ± —Et x | 3> \Et± —Et x \, the T± -B- T x 
transitions are the most important (see figure l2bj figure l2cl and figured]). The absolute 
energy difference between these levels is given by the total effective Zeeman splitting 
(see (ED): 

A to t= (g^BB cxt + ^-p). (37) 



2 

As the HF coupling intensity increases, so does the energy difference between the initial 
and the final states. Therefore, absorption processes become weaker with respect to the 
emission processes as At increases (see (EE}), allowing the nuclei spins to become more 
polarized (figure [3]). 

Figure 0] shows the energy levels of the DQD versus the external magnetic field. 
Electron- nuclei spin feedback is taken into account, and the external field is swept 
forwards and backwards. Figure Hal shows the energy levels for the smallest HF coupling 
considered (A L = 70 fieV). In this case: i) the induced Overhauser field is always parallel 
to the external field (see ([6])); and ii) the T±T X crossing occurs at B ext = B T t = 0T. 
Figure lib] shows that for the largest value of the HF coupling (At, — 90 fieV) considered 
here, the effect of including the feedback is to renormalize the electronic energy levels in 
such a way that the T±T X crossing occurs at larger absolute magnetic field than in the 
previous case. When sweeping forwards (backwards) from B ext < (B ext > 0) through 
-Bext = 0, there is a negative (positive) nuclei spin polarization built up (figure 13c]) . hence, 
it is still necessary to increase (decrease) the external field in order to compensate the 
accumulated Overhauser field and reach the T±T X crossing. This crossing occurs now 
for Btt > (Btt < 0) when sweeping forwards (backwards). Thus, figure [3b] shows 
a small region for positive values of the external field where the Overhauser field and 
the external field are antiparallel. Finally, figure I3c1 and figure 0b] show that, precisely 
in this region where the external and the induced fields are antiparallel, hysteresis is 
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observed as the external field is swept backwards. Notice that the hysteresis is observed 
only for the largest value of the HF coupling intensity considered. Moreover, the size of 
the hysteresis loop increases with the HF coupling intensity (not shown in the figures). 
Summarizing: i) The amount of polarization induced in the nuclei spins depends on the 
competition between absorption and emission spin-flip processes; ii) the nuclear spin 
polarization increases with the HF coupling; and iii) for large values of the HF coupling, 
the nuclear spin polarization versus the external magnetic field presents a bistable region 
where hysteresis is observed. In this region the external field and the Overhauser field 
are antiparallel. 

We will now analyze the leakage current behaviour versus the external magnetic 
field. Figure [3d] (without feedback), figure [3e] (with feedback) and figure [3f| (sweeping 
forwards and backwards with feedback) shows different behaviours for the current 
depending on the intensity of the HF coupling. The smallest HF coupling intensity 
considered (A L = 70 /ieV) presents a current dip around the T±T X crossing (figure H]), 
whereas the largest one (Al = 90 /xeV) presents a current peak. This behaviour can 
be understood recalling the current behaviour studied in the scheme given in f l2T]) . 
Briefly, this scheme defined a low and a high current regime, comparing the T x — > T± 
spin-flip rates (Wt±,tx) with the T x — > |±) tunnelling rates (T±^ Tx ). These rates are 
shown in figure |5j Figure [5a] (A L = 70 /ieV) shows that around the T±T X crossing 
Wt ± ,t x 3> T±,t x (w+ ^> (3 in diagram (12TT)). Therefore, it is more probable for electrons 
to spin-flip from T x — > T± than to tunnel from T x — Y |±) through the contact barrier, 
so a current dip is observed (low current regime). By contrast, figure [5c] {Al = 90 /zeV) 
shows that around the T±T X crossing, Wt ± ,t x < ^±,t x (w+ < /? in diagram ( )2T|) ). 
Therefore, it is more probable for electrons to tunnel from T x — > |±) than to spin- 
flip from T x — y T±, so a current is enhanced and a peak is observed (high current 
regime). Furthermore, Figure [3a] shows that for the smallest intensity considered for 
the HF coupling (Al = 70/ieV), the nuclear spin polarization goes to zero around the 
T±T X crossing, and so does T± jTx (see ( 131)) and figure [5aJ, thus, a current dip is observed 
(figure l3dl) . However, for the largest intensity considered (Al = 90/ieV), the nuclear 
spin polarization is finite around the T±T X crossing and T± )Tx > Wt ± ,tx (figure [5c]) , thus, 
a current peak is observed (figure l3dl) . 

Figure [3d] figure [3e] and figure Hf] show how the leakage current increases as the HF 
coupling intensity increases when the external field is close to zero. When increasing the 
HF coupling, the tunnelling rate T± ) t x (see (13^1) ) also increases (figure |5]) and, therefore, 
so does the current. For Al = 70 fieV and Al = 80 /ieV, a current drops to zero at 
B ext = Btt = 0, whereas for Al = 90 /ieV a current peak (figure l3d]) is observed. This 
peak occurs at B ext = B TT ^ when feedback is considered (figure |3e]). Moreover, 
figure |3e] shows that the effect of the electron nuclei spin feedback is only appreciable for 
the largest value of the HF coupling (Al = 90 /zeV). Therefore, only in this case current 
hysteresis is observed when sweeping backwards the external field (figure |3fj) . 
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Figure 3: a), b), c) Induced nuclear spin polarization versus external magnetic field, 
d), e), f) Leakage current versus external magnetic field. Parameters: T = 120 mK, 
t LR = 50/xeV, T L = T R = k B T/3, A R = 0.8A L , 7 = 5/ieV, N = 5 x 10 4 . Initial 
conditions: pi R (0) = 1 and -P(O) =0. A L = 70 yueV (solid, red online), A L = 80 /ieV 
(dashed, green online) and Al = 90 /xeV (dotted, blue online). Hysteresis is only 
observed for the largest value of the HF coupling (A L = 90/ieV). The leakage current 
shows a dip (peak) for the smallest (largest) value of the HF coupling. 




38 L , , 1 . , J t , ■ . 1 . . lad 
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(a) Al = 70 ^eV (b) A L = 90 ^eV 

Figure 4: Energy levels versus external magnetic field taking the spin electron-nuclei 
feedback into account and sweeping forwards and backwards B ext , for A L = 70 fieV (a) 
and A L = 90 /zeV (b). T + (solid, red online), T_ (dotted, blue online), T x (dashed, 
green online). Same parameters and initial conditions than in figure EJ The magnetic 
field range considered just includes the T±T X crossing. Singlet states are far away in 
energy. As in figure [3l hysteresis is only observed for the largest value of the HF coupling 
(A L = 90/ieV). 
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(a) A L = 70 /jteV. (b) A L = 80 ^eV. (c) A L = 90 ueV. 

Figure 5: Spin-flip rates (a;*, see f lT9|) ; solid, red online) and tunnelling rates through 
contact barriers (/3, see f lT9|) ; dashed, blue online) versus external magnetic field. Same 
parameters and initial conditions than in figure [3J In figure |Sa] and figure [5b] there 
are two regions: i) (3 > u;+, namely, electron tunnelling from T x to the right lead is 
more effective than electron spin- flip from T x to T± triplets (high current regime); and 
ii) j3 < a;", namely, electron tunnelling from T x to the right lead is less effective than 
electron spin-flip from T x to T± triplets (low current regime). In this case, the current 
shows a dip around zero external magnetic field (figure [3]). In figure l5c| f3 > cu+ always. 
In this case, the current shows a peak (figure [3]). 



3.1.2. Simplified model around the triplet-triplet crossing In the regime described 
in section 13.1. 1[ the relations \E T± — E s+ \ 3> 7 and \E T± — Eg_\ 3> 7 are well satisfied, 
so the rates u± that appear in the reduced rate equations (1201) can be safely neglected 
(see dl9]) and ([29}) and ([20]) becomes: 

P = co^ + co^p Tx - ujKp (38) 

Considering the stationary limit, we obtain the following equations for the stationary 
solutions: 

4+ -4 

— 2 Vt (39b) 



3 + 

where p t = co^/co^. The rates and (3 are complicated functions of P (see (TT9|) . 
( |2^|) and (I30al) ). thus, ( 139b[) cannot be solved analytically We have obtained numerical 
solutions for £> ext = (figure [6]) case, although f|39ap and f[39b[) hold also for finite 
external fields. Figure [6] shows that increasing the HF coupling intensity produces a 
bifurcation on the induced nuclei spin polarization. For A L < 80.43 /ieV, the nuclei spins 
have one stable solution at P = 0, namely, the nuclei spins are fully depolarized (solid 
line). However, for Al > 80.43 /ieV, the P = solution becomes unstable (dashed line), 
and two stable solutions (solid lines), with the same absolute value but with opposite 
signs appear. This behaviour is the same found in the full numerical solution (figure l3al) . 
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We have found that the system undergoes the bifurcation as the slope of the linear term 
of the expansion around P = of the right hand side of f!39bj) is varied [39J. The slope 
is given by: 

s= If^- 4 V (40) 



3 \k B T 

When s < 1, P — is the only stable solution, while when s > 1 the system presents the 
two stable solutions mentioned above. The bifurcation occurs at s = 1. Putting into (140]) 
the parameters for which figure [6] is obtained we find that the bifurcation takes place 
when Al = 80.43 fieV. Therefore, ( [40]) provides an expression that relates the hyperfine 
coupling intensity with the temperature and the bifurcation which, in addition, can be 
observed through the hysteresis plots measured for the current through the DQD versus 
the external field [3j|lJ[6j[T0]. Finally, we find that the agreement between the induced 
nuclei spin polarization obtained with the full calculation for S ext = 0, and the one 
obtained with this simplified model is very good. 

From the current through the right contact barrier Ir = (T +) t x + T-,t x )Pt x 
(see (135p ). and from (J39aj) we obtain the following expression for the current through 
the DQD: 

Notice that (3 = when P = (see ffl9|) ). thus, the current will be zero when the nuclei 
spins are fully depolarized. In the previous discussion, we have shown that the induced 
nuclei spin polarization shows a bifurcation when increasing the HF coupling intensity. 
Therefore, in the range of Al where the nuclei spins are fully depolarized (figure [6]) no 
current flows through the DQD at B ext = 0. However, in the range of Al where the 
nuclei have a non zero spin polarization (figure [6]), a finite current flows through the 
DQD at -Bext = 0, since I(B cxt = 0) oc P 2 to lowest order in P. Furthermore, ( T4T1) shows 
that the current depends on the ratio These rates are shown in figure [5j and 

motivate the physical picture we have used to understand the transition from a current 
dip to a current peak. Therefore, this simplified model shows that the bifurcation 
obtained for the induced nuclei spin polarization together with the ratio /3/ w +> describe 
the transition from a current dip to a current peak (figure l3d|) . 

3.1.3. Large magnetic fields In this section, we consider larger magnetic fields than 
in the previous case in order to account for the ST crossings. This means sweeping 
the external field from B cxt < -Bs ± ,t t to B ext > Bs ± ,t ± (see schematic figure [2a| 
and figure l2dl) . Therefore all three level crossings (S±T T , S±T± and T±T X ) will be 
considered (see figure ITal and figure l7bl) . Figure [7] shows the DQD energy levels, the 
induced nuclear spin polarization and the leakage current through the DQD versus the 
external magnetic field, for Al = 70 fieV and Al = 90 /ieV. In this case all figures 
include feedback and the external field is swept forwards and backwards. 
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Figure 6: Solutions for ( I39bl) for B ext = 0. The induced nuclei spin polarization 
presents a bifurcation. When Al < 80.43 /xeV, there is only one stable solution, which 
corresponds to having the nuclei spins fully depolarized (solid line). However, when 
A L > 80.43 /ieV, P = becomes an unstable solution (dashed line) and two stable 
solutions show up (solid lines). Same parameters as in figure [31 

Nuclear spin polarization versus the external field is shown in figure [7c] (A^ = 
70 /ieV) and figure [7d] (A L = 90/zeV). The external field sweeping starts for B cxt < 
Bs ± ,t 1 :- Initially, the nuclei spins are completely depolarized. The DQD energy level 
distribution shows that T + triplet is the ground state (see figure [7a] and figure l7bj) . In 
this case, all spin-flip emission rates polarize negatively the nuclei spins (see (136ajl ) and, 
hence, they become almost completely negatively polarized (P ~ —1). Furthermore, 
since the probability of finding a nuclei with spin +1/2 is nearly zero when P ~ — 1 
(see ( |24l) ). all emission rates become approximately zero (see (1291)). and the nuclei 
spin polarization remains roughly constant until _B ex t reaches the S±T T crossings. 
Only close enough to the S±T T level crossings, spin-flip absorption processes become 
significant and the nuclei spins dramatically depolarize (figure [Tel and figure l7dl) . In the 
Bs±,t 1 : < B cy± < Bs ± ,t ± range of magnetic fields (schematic figure l2bl and figure l2cl) . 
several spin-flip rates compete in order to polarize the nuclei spins in opposite directions 
(see (136bj) and (I36cj) ) . When B s± ^ < B cxt < B TT (scheme in figure l2dl) . most of the 
spin-flip emission rates polarize negatively the nuclei spins (see (I36bl) ). and the resulting 
nuclear spin polarization is negative in this region. However, as the levels approach the 
T±T X crossings, spin-flip absorption processes become more relevant, and the nuclei 
spins slowly depolarize. When B TT < B cxt < B s±) t± (scheme in figure [2c]) . most of the 
spin-flip emission rates polarize positively the nuclei spins (see (I36cl) ). and the nuclear 
spin polarization is positive in this region. Finally, when B ey± > Bs ± ,t ± , the T_ triplet 
is the ground state (scheme in figure l2dl) . and (136d|) shows that all spin-flip emission 
rates polarize positively the nuclei spins, thus, they become almost completely positively 
polarized (P ~ 1). The probability of finding a nuclei with spin —1/2 is nearly zero as 
P ~ 1 (see ( 1241) ). thus, all emission rates become approximately zero (see ([29]) ). and the 
nuclei spin polarization remains roughly constant. 

Figure [7c] and figure [7d] show that the feedback between electron and nuclei spins 
produces hysteresis in the nuclear polarization around the ST-crossings. Recall that for 
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small magnetic fields (section [3.1. ip . where only T±T X crossings participate, hysteresis 
was not observed for small HF intensities, but only the largest HF coupling intensity 
considered here (Al = 90/zeV). For larger external magnetic fields, however, hysteresis 
shows up at ST crossings even for the smallest HF intensity considered (A L = 70 fieV). 
Finally, as for small magnetic fields, hysteresis is larger as the HF coupling increases. 

The leakage current through the DQD versus the external field is shown in figure [Tel 
(Al — 70 /xeV) and figure [7fl (A L = 90/xeV). The current presents three peaks, each of 
them corresponding to one of the three possible level crossings. When B cxt < B$ ± t t 
(B ext > B s± t±), the current is zero because electrons are trapped in the T + (T_) triplet 
state (figure ITal and figure l7bl . thus, in these ranges of the external field spin-flip emission 
rates are nearly zero. In between these crossings (when B s± ^ < B ext < B s±jT± ), current 
is strongly quenched but nevertheless finite. This case has been discussed in section 13.1.1 1 
Finally, the feedback between the nuclei and the electron spins also produces hysteresis 
in the current around the ST crossings, when sweeping the external field forwards and 
backwards. 

3.2. Dependence on the interdot tunnelling strength 

In this section, we show the leakage current and the induced nuclear spin polarization 
dependence on the interdot tunnelling intensity. Interdot tunnel varies a lot from one 
experiment to another and can be externally tuned. For instance, in [2] it is estimated to 
be around 30 /xeV, while in [3] around 0.2 /xeV. This justifies the large difference between 
the two values that we have chosen. We consider small external magnetic fields, and 
two different interdot tunnel values (figure |H_a] and figure l8bl) . For the largest value of 
the interdot tunnel (tLR = 50 /xeV) only the T±T X crossing participate in the current. 
However, for the smallest value of the interdot tunnel (t LR = 0.01 /xeV), all crossings: 
S±,T T , T±T X and S±,T± participate in the current, for the same range of the external 
field. Figure [8] shows also the induced nuclear spin polarization (figure IHcl and figure l8d|) 
and the leakage current through the DQD (figure [He] and figure (HJ) versus the external 
magnetic field for the two interdot tunnelling intensities chosen: tLR = 50 /xeV and 
t LR = 0.01 fieV, and for A L = 80/xeV. 

Figure [H_c] and figure [8d] show that for the smallest interdot tunnelling (tLR = 
0.01 /ieV), the behaviour of the polarization versus the external magnetic field is 
smoother than for the largest interdot tunnelling considered (tLR = 50/xeV), due to 
the stronger competition between the energy absorption and emission processes in the 
former case. For the small interdot tunnelling, absorption is more efficient than for the 
large interdot tunnelling value. 

Figure [He] and figure [8f] show that decreasing the interdot tunnelling, the current 
versus the magnetic field presents again a transition from a dip to a peak. As in the 
previous case, where we discussed the current behaviour as a function of the intensity of 
the HF interaction, a dip or a peak features observed in the current can be understood 
comparing the spin-flip rate between T± and T x and the tunnelling rate through the 
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Figure 7: a), b) Energy levels versus external magnetic field. T + (solid, red online), 
T_ (dotted, blue online), T x , S+, S- (dashed, green online); c), d) induced nuclear spin 
polarization, and e), f) leakage current versus external magnetic field sweeping forwards 
(solid, red online) and backwards (dashed, green online). Same parameters and initial 
conditions than in figure [3l The range of magnetic fields considered includes both ST and 
the T±T X crossings. Hysteresis is observed for both values of the HF coupling intensity 
at the ST crossings. Recall that in figure [3l where only the T± crossing was shown, 
hysteresis was only found for the largest value of the HF intensity (A^ = 90 /xeV). The 
current shows now three peaks corresponding mainly to each of the levels crossings. 



barrier contact between T x and |±) states (figure l5bl) . However, unlike in the previous 
case, where the HF coupling had to be increased to observe the transition, now, in order 
to go from a dip to a peak, the interdot tunnelling must decrease. 
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Figure 8: a), b) Energy levels versus external magnetic field. T + (solid, red online), 
T_ (dotted, blue online), T x , S+, S- (dashed, green online); c), d) induced nuclear spin 
polarization, and e), f) leakage current versus external magnetic field. Al = 80 fieV. 
The remaining parameters and the initial conditions are the same as in figure [3j The 
current shows a dip (peak) for the largest (smallest) value of the interdot tunnel intensity. 
Although feedback between the electron and the nuclear spin polarization is taken into 
account in these figures, hysteresis is not observed for this set of parameters. 



4. Conclusions 

We have studied the leakage current through a coherently coupled DQD in SB regime. 
Spin relaxation due to HF interaction between the spins of the electrons in the DQD 
and the nuclei spins lifts SB producing leakage current. Moreover, the spin interaction 
between electrons tunnelling through the DQD and nuclei, induces dynamical nuclear 
spin polarization that is in general non-negligible. We have investigated the behaviour 
of both the induced nuclear spin polarization and the leakage current as a function 
of an external magnetic field. Our three main results are: i) the leakage current 
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shows a dip or a peak depending on the intensities of both the HF interaction and the 
interdot tunnel strength. We have shown that for large (small) HF coupling (interdot 
tunnelling) intensities the current shows a peak. On the contrary, for small (large) 
HF coupling (interdot tunnelling) intensities the current shows a dip. Large (small) 
HF couplings (interdot tunnelling) indicate strong mixing between the singlet-triplet 
subspaces, namely, the effective Zeeman splitting difference between the dots is non- 
negligible with respect to the exchange energy. In this case the leakage current shows 
a peak. On the contrary, small (large) HF couplings (interdot tunnelling) indicate 
a weak mixing between the singlet-triplet subspaces, the Zeeman splitting difference 
between dots becomes negligible with respect to the exchange energy, and the system 
is mostly blocked in the triplet subspace. In this case the leakage current shows a dip. 
The crossover from a dip to a peak is, thus, obtained by increasing (decreasing) the HF 
interaction (interdot tunnelling strength), ii) For a wide external magnetic field sweeping 
range, we have shown that the leakage current shows three main peaks. Two satellite 
peaks corresponding to the ST crossings and a central peak corresponding to the triplets 
crossing, iii) We have observed hysteresis in both the leakage current and the induced 
nuclear spin polarization as a function of an external magnetic field. This hysteretic 
behaviour is a consequence of the dynamical nuclear spin polarization interacting with 
the electron spin that tunnel through the DQD structure. Finally, we have shown that 
the size of the hysteresis region strongly depends on the HF interaction intensity and 
the interdot tunnel strength. 

Our results are a contribution to ongoing efforts at understanding and controlling 
spin relaxation in qubits, which limitate their performance for quantum information and 
quantum computation purposes. 
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